A possible experimental test of the thermodynamic approach to granular media 
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We study the steady state distribution of the energy of the Sherrington-Kirkpatrick model driven 
by a tapping mechanism which mimics the mechanically driven dynamics of granular media. The 
dynamics consists of two phases: a zero temperature relaxation phase which leads the system to 
a metastable state, then a tapping which excites the system thus reactivating the relaxational 
dynamics. Numerically we investigate whether the distribution of the energies of the blocked states 
obtained agrees with a simple canonical form of the Edwards measure. It is found that this canonical 
measure is in good agreement with the dynamically measured energy distribution. A possible 
experimental test of the Edwards measure based on the study here is proposed. 

PACS numbers; 05.20, 75.10 Nr, 81.05 Rm. 



Complex systems such as granular media possess a 
large number of metastable or blocked configurations. 
When a granular medium is shaken it quickly relaxes into 
a blocked configuration, a subsequent shake or tap will 
lead it to another blocked or jammed state and so on. 
If the driving mechanism is held constant one expects 
the system to enter into a quasi-equilibrium stationary 
state. Various driving mechanisms can be investigated 
experimentally, such as vertical tapping Q and horizon- 
tal shaking . In granular media and other complex sys- 
tems such as spin glasses, the entropy of these blocked 
states is extensive in the system size, and hence it has 
been proposed that one may use a thermodynamic mea- 
sure over blocked states to describe this steady state. The 
simplest proposition is that the system is characterized 
by a number of quantities which are fixed on average and 
then the measure on the steady state is obtained from the 
maximum entropy state (on blocked states) with the rel- 
evant macroscopic quantities fixed This simple idea 
has recently been investigated in a wide range of systems 
and has been shown to be relatively successful. Various 
tests of the applicability of these thermodynamic ideas 
have been carried out and although some confirmation 
has been made in more realistic sheared granular systems 
Q , most work has been carried out on simpler model sys- 
tems which one hopes capture the basic physics of gran- 
ular media. The Edwards flat measure has been shown 
to be of predictive value in some simple one dimensional 
lattice models [^|| , but there are clearly examples where 
the approach fails 0. However, recently it was shown 
that more sophisticated versions of the Edwards measure 
introducing ensembles with several quantities fixed on av- 
erage can remedy the deficiencies of the basic measure in 
these cases Other toy models that have been ana- 
lyzed are lattice based models with kinetic constraints in 
higher dimensions ||^,|l^ and also spin glass models pl| ] 
where non-thermal driving is used to move the system 
between blocked states. 

Even if it is not expected to be exact, many sys- 



tems may be described to a good engineering level by 
these measures. Given the difficulty of the analysis of 
the highly nonlocal dynamics in these systems this is an 
important step toward understanding their steady state 
regimes. There is no clear ergodicity in these systems 
and no detailed balance as in usual statistical mechan- 
ics. Edwards argued that a system might conceivably 
explore blocked configurations in a flat manner if the 
driving involved extensive manipulations, meaning the 
displacement of a macroscopic number of particles, for 
example shaking, stirring or pouring granular media. An 
interesting consequence of the applicability of thermody- 
namic ideas is that one may describe phase transitions 
in these driven systems [|l2|. However, even considering 
the success of the Edwards measure in describing various 
simple models, evidence in realistic granular media is still 
lacking. In this Letter, we investigate whether, at fixed 
tapping rate (to be defined later) , the states explored dy- 
namically obey a form of Boltzmann distribution. The 
results presented here are quite striking, despite a lack of 
detailed balance we shall see that a Boltzmann distribu- 
tion excellently describes the histogram of the energies 
of the blocked states visited during the tapping. Moti- 
vated by these results, we propose a generic and simple 
experimental test of the Edwards measure, which should 
be feasible in a wide range of driven granular systems. 

As mentioned above, a good theoretical and numer- 
ical testing ground for this thermodynamic approach 
to granular media are spin glasses. Spin glasses have 
an exponentially large (with the system size) number 
of metastable states and hence an extensive entropy of 
blocked states as do granular media. The definition of a 
blocked state in a spin glass simulated on a computer de- 
pends of course on the local dynamics. Under single spin 
flip dynamics a metastable state is one where flipping 
any single spin increases the energy, it is thus a blocked 
state under any single spin flip Monte Carlo dynamics. 
Various spin glass models have been studied to explore 
the accuracy of the Edwards measure as a function of 
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the relaxational dynamics and the tapping mechanism 
p| . Here we shall explore the driven dynamics of the 
Sherrington-Kirkpatrick (SK) spin glass model de- 
fined by the Hamiltonian H — ^I^kj Jij^i^j^ where 
the Si {I < i < N) are Ising spins and the interaction 
Jij are independent Gaussian random variables of zero 
mean and variance J?- = 1/N. 

The tapping mechanism introduced in spin glass mod- 
els involves two steps: 1) Relaxation: the system evolves 
under zero temperature Glauber dynamics until it is 
blocked in a metastable state; 2) The system is tapped 
with tapping parameter p, that is to say each spin is 
flipped in parallel with probability p - this is the exten- 
sive manipulation as at each tap on average pN spins are 
flipped. Given a tapping rate p, it is observed that the 
system reaches a steady state regime with a fixed aver- 
age energy per spin e (the energy being measured at the 
end of step 1) above). If we postulate the canonical form 
of the Edwards measure for the quasi-equilibrium distri- 
bution of the driven dynamics (observed on the blocked 
states) one has the steady state distribution of energy 
given by 
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where Nms{E) is the number of metastable states of en- 
ergy E, (3 is the inverse Edwards temperature, which can 
be thought of as a Lagrange multiplier fixing the aver- 
age energy per spin, and Z = 7Vms(-E) exp(— /3i?) 
is the canonical partition function. We note that a 
fully flat measure with fixed energy per spin would be 
Pa = exp{—/3Ea)/Z, where Ea is the energy of the 
blocked state a. In this paper we shall investigate only 
the weaker form of the measure Eq. (^). Such an in- 
vestigation has been carried out in a model of particle 
deposition , but the originality of the approach of this 
Letter is that one can partially test Eq. (|^) without 
knowing the Edwards entropy. This is crucial, since it 
opens up the possibility of looking at a wider range of 
granular models and even experiments. 

The hypothesis Eq. (^ will be tested in the follow- 
ing manner. We fix a bin size for the energies AE and 
by exact enumeration we compute N]^Js{E) for a given 
system and realization of the disorder. The exact enu- 
meration means that we can explore system sizes up to 
N = 30. Then the tapping dynamics is applied to the 
same system and the dynamical histogram N]j{E,p) of 
the energy with the same bin size AE is computed over 
10^ taps after the system has been tapped 10^ times to 
reach the steady state. If the Edwards measure in the 
form Eq. (0) is valid then one should obtain 
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where C is a constant which should be independent of E. 
If Eq. (P) is valid then an interesting question to ask is 



what is the dependence of (3 on a) the tapping strength 
p, b) the size of the system N c) the realization of the 
disorder. It is interesting to see if the tapping strength 
p determines a unique temperature in the limit of large 
system sizes, for systems with the same statistical dis- 
tribution of couplings i.e. is /3 self averaging? We shall 
see that the results obtained here for small system sizes 
suggest that Eq. (||) is obeyed on increasing the sys- 
tem size and that the parameter [3 becomes independent 
of the system size and the disorder realization for large 
systems. 

In order to see what happens for larger systems, where 
the enumeration of metastable states is no longer feasi- 
ble and where the thermodynamic approach should be 
expected to work better, one may compare the dynami- 
cal distribution N]:i{E,p) measured for the steady states 
at different tapping rates. If the Edwards measure holds, 
i.e. N]j{E,p) cx pedw{E), then for any two tapping 
strengths p and p' we should find that 
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or equivalently 



exp {-{(3{p) - f3{p'))E + f{N,p,p')) (3) 



HNd{E,p)) - \n{ND{E,p')) 
N 

-{(3{p)-f3{p'j)e+^f{N,p,p') (4) 



where e — E/N is the energy per spin and f{N,p,p') is 
independent of E. This constitutes a weak test of the Ed- 
wards hypothesis but has the merit of being practicable 
for large system sizes. 

A final test is to use the analytic, but annealed calcula- 
tion, of the average number of metastable states at fixed 
energy, Nms (e) (the overline indicates the disorder aver- 
age over realizations) , which gives in the thermodynamic 
limit NMs{e) « eMNSlisie)), where 0^ 
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is the annealed entropy per spin of metastable states of 
energy e per spin. If the Edwards measure holds then we 
should find that 
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= -I3{p)e + C 



(6) 



We further expect that the dynamical entropy per spin 
Sd{^,p) = ln(iV£i(e,p))/A^ and the entropy of metastable 
states per spin S]\is{e) = ln(iVj\,fs(e,p))/A^ are self av- 
eraging and thus Eq. ( | 6[) m ay be written for large N 
as r(e) = Soi^.p) ~ Susie) = -/?(p )£ + C. Unfortu- 
nately no exact calculation for S'Ms(e) exists, however it 
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is known that (e) = Snsie) for > -0.672 [||. 
Furthermore this annealed approximation seems to be 
good, though not exact, to even lower energies as con- 
firmed by a replica symmetric calculation [|l6|. Hence at 
energies above Ec we have the formula 

r,{e)^S^{^~S-,js{e)^-m^ + C (7) 

which should also work well approximately at energies 
below, but not too far from, ec- We remark here that by 
Jensen's inequality 5^/5 (e) > S'Ms(e)- 
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FIG. 1. The computed values of r{E) for the randomized 
spin fiip relaxation at p = 0.3 with a bin size of 0.025A'^ for 
systems sizes A'' of 15 (circles), 20 (squares), 25 (diamonds) 
and 30 (triangles). The straight lines are linear fits to guide 
the eye. 

We now discuss the results of these numerical tests. 
Systems between size 15 and 30 were analyzed. The 
numerically measured value of r{E) was computed from 
the static metastable state energy histogram and the his- 
togram obtained by tapping the system 10^ times. This 
was done for values of p, 0.1, 0.2, 0.3, 0.4 and 0.5. In- 
creasing the number of taps did not change the dynamical 
histogram so we can be sure that we are measuring the 
steady state (this is natural given the small system sizes). 
It is possible that with very low tapping rates the time 
to reach the steady state even for a system of 30 spins 
becomes large. As an example the results for p = 0.3 are 
shown in Fig. (|^). We sec that for small systems sizes 
the points of the numerically measured r(E) are scat- 
tered about the straight line fit predicted by Eq. (|^). 
However as the system size is increased the straight line 
fit appears excellent. The deviations from the straight 
line fit also appear to be non systematic. Similarly good 
fits are also obtained for all the tapping rates studied. 
In the case of the SK model, a Boltzmann distribution 
of the energies of the metastable states visited in the 



steady state regime appears to depend simply on hav- 
ing a large enough system size. It is generally believed 
that the Edwards measure should be more accurate in 
gently driven systems, here the weak form characterized 
by Eq. (||) seems equally valid at all tapping rates. We 
also note that the slope of the straight lines appear to be 
saturating to a constant value of /3(p) on increasing TV, 
suggesting that in the thermodynamic limit the Edwards 
temperature is fully characterized by the tapping rate p. 

From this first numerical study one may conclude that 
if the weak form of the Edwards measure holds then it 
relies on the system being sufficiently large. The re- 
sults above however seem encouraging, especially when 
one takes into account that the average total number of 
metastable states in a system of size N exp(0.1992Af) 
|]l^ , which is only about 400 for a system of size 30. 

We now carry out simulations for a system of size 200, 
here the average number of metastable states is about 
2 X 10^"^. The bin size is also refined to 0.00257V. Naively 
one would expect the thermodynamics approach to work 
better in this limit as the phase space of metastable states 
is much bigger. Also there are many more states in each 
energy bin and hence it is really the weak form of the 
Edwards measure that is tested. Shown on Fig. (||) is 
the numerically computed c{p,p,' e) from a simulation of 
the same system tapped at values of p, 0.1, 0.2, 0.3, 0.4 
and 0.5. The system was first equilibrated by tapping 
5000 times and the dynamical histogram was constructed 
over 5 x 10^ taps. One sees from Fig. (2) that the val- 
ues of c{p,p',e) are to a good degree of accuracy on a 
straight line, as predicted by Eq. (^). The deviations at 
low energies, where the sampling is lowest, appear non- 
systematic. 
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FIG. 2. The numerically computed values of c{p,p'e) from 
simulations of a system of size A'^ = 200 spins. The curves 
have been vertically shifted for clarity. Shown in descending 
order are the results for {p,p'): (0.2,0.1), (0.3,0.2), (0.4,0.3), 
(0.5,0.4). 
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The energy histograms generated by the precedent 
simulations were used to compute ^^(e), the results are 
shown in Fig. (|^). The curves have been shifted verti- 
cally for clarity. A straight line fit was performed in the 
region e > Ec and we see that in this region the fit is ex- 
cellent. In addition the values obtained from the fits for 
f3{p) are compatible with those obtained for the systems 
of size = 30. The lower energy part of the curves also 
appear linear down to around e = —0.72. The replica 
symmetric calculation of also visibly (on the same 
scale) departs from the annealed calculation near this 
energy, although we emphasize that this is not an exact 
calculation but can be expected to be closer to the real 
result than the annealed one. The lower tapping rates 
explore the lower energy regime while the higher rates 
explore the higher energy regime. The deviations above 

from the annealed entropy are at the high energy end 
of the numerical histogram where the sampling is small- 
est. 
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FIG. 3. The calculation of r-a(e) for a system with A'^ = 200 
and data taken over 5 x 10® taps. Shown in ascending order 
are the plots for p: 0.1, 0.2, 0.3, 0.4 and 0.5. The straight line 
fits are shown in the region e > Ec (to the right of the vertical 
line) . 

To summarize we see that the assumption of an effec- 
tive Boltzmann distribution for the energy of the states 
explored by tapping dynamics describes extremely well 
the numerically obtained results. This description seems 
to work better on increasing the system size for small sys- 
tems and seems to be valid for larger systems via slightly 
more indirect tests. Why it works is rather mysterious 
and if it is exact it would be an extremely useful method 
to numerically map out the distribution of metastable 
states in various spin glass models. The relaxational 
dynamics used was random update which uses a lot of 
computer time for these simulations. Sequential update 
is much quicker for these simulations and we have car- 
ried out such simulations and found to within very small 



deviations the same results. If one could prove the ap- 
plicability of the Edwards measure for these types of dy- 
namics one would have an extremely powerful method to 
explore metastable states and inherent states in glassy 
systems. The standard method employing an auxiliary 
Hamiltonian is much more time consuming - although 
classical statistical mechanics tells us that it will give the 
right result. 

In the original proposition of Edwards in the context 
of granular media, the steady state volume fraction V 
occupied by a driven granular media is described in the 
canonical approach by 



PDYN{V)^exp{~^ + S{V)) 



(8) 



where X is the compactivity and S{V) is the entropy of 
blocked states occupying volume V. The problem of cal- 
culating S{V) in a real system seems formidable. How- 
ever it should be possible to test the prediction of Eq. 
(^) experimentally by comparing the histograms of V 
for systems driven at different driving amplitudes. This 
would be a crucial test, though we emphasize again, not 
a demonstration of the validity of the Edwards measure. 
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